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Abstract 



We examine the validity of the hypothesis of self-similarity in systems 
coarsening under the driving force of interface energy reduction in which three 
dimensional particles are intersected by a one or two dimensional diffusion 
matrix. In both cases, solute fluxes onto the surface of the particles, assumed 
spherical, depend on both particle radius and inter-particle distance. We 
argue that overall mass conservation requires independent scalings for particle 
sizes and inter-particle distances under magnification of the structure, and 
predict power law growth for the average particle size in the case of a one 
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dimensional matrix (3D/1D), and a weak breakdown of self-similarity in the 
two dimensional case (3D/2D). Numerical calculations confirm our predictions 
regarding self-similarity and power law growth of average particle size with 
an exponent 1/7 for the 3D/1D case, and provide evidence for the existence 
of logarithmic factors in the laws of boundary motion for the 3D/2D case. 
The latter indicate a weak breakdown of self-similarity. 
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Typeset using REVT^ 



I. INTRODUCTION 



The purpose of this paper is to re-examine the vahdity of the hypothesis of self-similarity 
in coarsening systems in which particles of dimension Dp are intersected by a diffusion ma- 
trix of different spatial dimensionality D. The discussion is limited to coarsening driven by 
interface energy reduction, the particles are assumed to remain spherical (e.g., by surface 
diffusion), and particle migration is neglected. We focus on two specific cases: three dimen- 
sional particles intersected by either a two (3D/2D) or one (3D/1D) dimensional diffusion 
matrix respectively. 

The statistical self-similarity hypothesis with one scaling length asserts that after a possi- 
ble transient, consecutive configurations of the coarsening structure are geometrically similar 
in a statistical sense As a consequence, any parameter of the structure that is invariant 
under a uniform magnification is also independent of time. This hypothesis, together with 
the laws of boundary motion for a specific system (and their scaling under uniform magni- 
fication) are sufficient to obtain the equation of motion for any linear scale of the structure 
(e.g., the average particle radius {R{t)) for an ensemble of coarsening spherical domains). 

Self-similarity with a single scaling length as stated above is consistent with conservation 
of mass (volume) only in systems for which Dp = D. Thus conservation of particle mass 
requires. 



where is the number of particles per unit (general) area of substrate and dav = I/'^a 
may be regarded as an average spacing between particles. If self-similarity holds for the 
particle size distribution, then 




'av 



(1) 
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(^Ri^p^ = const. {Rf" 



(2) 



so that Eq. (|I]) can be written 



(R) 



Dp/D 



d, 



= const. 



(3) 



•av 



Self-similarity with one scaling length requires {R) /dav = const. But according to Eq. 
this condition holds only if Dp = D. Therefore, self-similarity with a single scaling length 
cannot hold when Dp ^ D. Unfortunately, self-similarity with one scaling length was used 
in a section of Ref. to deduce erroneously exponents of particle size in a 3D/2D and a 
3D/1D system (t^^'^ and t^^^ respectively). 

If Dp 7^ D, self-similarity with two scaling lengths (e.g., (R) and dav) is still possible 
and is consistent with Eq. (^. In this case, the self-similarity hypothesis would assert that 
consecutive configurations of the system are statistically equivalent to those obtained by 
magnifying the original system using two scale factors, one for particle sizes and a second 
for inter-particle spacings; the two scale factors are related through Eq. (j^) to conserve 
mass. We shall present evidence that this occurs in the 3D/1D case and leads to a particle 
growth exponent of 1/7. 

In the 3D/2D system, we will present numerical evidence that a necessary condition 
for self similarity is not met. This condition, established in Appendix A, shows that if 
the distribution of reduced particle sizes is scale invariant (which it must be in any self- 
similar regime), then the ratio of the expected growth rates of any two particles must be 
scale invariant under magnification. We show numerically that this condition is not met 
in the case of a configuration with a large number of interacting particles of unequal sizes, 
randomly distributed in space with a Laplacian concentration field satisfying self-consistently 
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determined mean field boundary conditions at large distances from the particles. 

Support for the existence of self-similarity is quite strong in a large variety of physical 
systems for which Dp = D (see, e.g., reviews in In the case of three dimensional 

particles growing by diffusion on a two dimensional substrate (the 3D/2D case), experiments 
and theoretical investigations have been recently reviewed in |^. Briefly, the system of 
coarsening particles reaches (at least approximately) a scale invariant state in which the 
average particle radius grows in time as {R{t)) oc t^^^ [§. Such a growth law has been 
also predicted on the basis of mean field analyses [§-0. We shall present evidence that 



self-similarity with this growth law cannot be strictly true. 

In Section we use the law of boundary motion in a coarsening system, together with 



self-similarity, where applicable, to discuss growth exponents. In Section |T|, we present 
numerical evidence that in the 3D/1D system, the particle radii and inter-particle distances 
scale separately to obey Eq. (|^) with {R) oc t^/'' and dav oc t^^'^ . Section |V| addresses the 
3D/2D case and presents numerical evidence for the existence of a logarithmic factor in 
the laws of boundary motion that involves ratios of particle sizes to inter-particle distances, 
and indirect evidence that, contrary to the classical case of coarsening of two dimensional 
particles on a two dimensional substrate, the logarithmic factor leads to a (weak) breakdown 
of self-similarity. This breakdown would manifest itself in the existence of effective coarsening 
exponents that change very slowly in time. 

II. SELF-SIMILARITY AND GROWTH EXPONENTS 



In this section we present the derivation of the growth exponent for the case Dp = 3,D = 
1 and discuss the complication that arises in the case of Dp = 3, D = 2. The discussion is 
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based on the following theorem: If the reduced particle size distribution is independent of 
scale, that is, if 

n{R,t) = ^P{R/{R)), (4) 

where n{R,t)dR is the number of particles per unit "area" in the system with (volume 
equivalent) radii in dR, (R) is the average particle size and P{x), in which x = R/ (R), is 
the reduced particle size distribution function, then it is shown in the Appendix that 

(m) = ^-^G{xl (5) 

where (R\R^ is the expected value of dR/ dt for particles of radius R. The expression for G{x) 
given in the appendix shows that as x increases from zero, G increases from negative values, 
corresponding to particles that shrink on the average, to positive values corresponding to 
particles that grow on the average; G vanishes at a value x = Xc, corresponding to particles 
that, on the average, do not change size. Eq. (|^) shows that the ratio of the expected 
value of dR/dt for particles of two different sizes is independent of scale as stated in the 
introduction. We will use this property (which we will refer to as the ratio test) to argue 
that the 3D/2D system cannot be strictly scale invariant based on the numerical evidence 



presented in Section IV 



If self-similarity holds (with one or more scaling lengths), the growth exponent may be 
determined from the scaling of d {R) /dt. But Eq. (^ shows that d {R) /dt scales as (^R\R 
which in turn scales as dR/dt for a single particle. Hence the growth exponent may be 
determined from the scaling of dR/dt. 

To discuss the scaling of dR/dt, we consider for simplicity an ensemble of precipitate 
particles embedded in a matrix, such that their growth or dissolution is limited by diffusion 



of one of the species (solute) through the matrix (the discussion can be easily generalized 
to other systems in which coarsening is driven by surface free energy reduction [00). Let 
c be the concentration of solute in the matrix, and assume that c in the matrix is much 
smaller than the (constant) concentration in the precipitate phase Cp, so that Cp — c can 
be approximated by Cp [|T^. Conservation of solute mass requires that for each precipitate 
particle, 

where j„ is the normal component of the solute flux along the outward normal to the 
surface, is the surface of a particle of radius R in Dp dimensions and is, for example, 
Sd^ = Ti^^/'^R^p^^ /V{Dp/2) for a hemisphere, and the integral is taken over the surface of 
the precipitate particle. 

For the 3D/1D case, the collection area over which is nonzero is of microscopic size 
and independent of scale. Let all particle radii be multiplied by A. Then the flux j„ scales as 
a gradient which scales as l/{\dav)-, the factor A arising from the Gibbs- Thomson equation. 
But, from Eq. (|^) X^/dav = const. Hence Eq. (^) shows that <^-R|-R^ scales as A~®. Therefore, 
if self-similarity holds, Eq. (|^) shows that d (R) / dt scales as A~®, or as (R)^^ so that 

^ = const. (7) 

which integrates to {R) oc t^/'' and dav oc t^/^. Evidence supporting this result in presented 



in Section III 



To discuss the case of a two dimensional matrix we first present an expression for dR/dt 
in the spirit of a mean field treatment. Consider a disk of radius R located at the origin, 
and assume quasi-steady diffusion in the matrix (V^c = 0), with boundary conditions c(r = 



^) = at some cut-off distance away from the center, and 

c(r = R)=co{l + ^y (8) 

at the disk's boundary, where Cq is the solute concentration in the matrix at coexistence, and 
r is the capillary length. Then the rate of change of volume of the particle is proportional 
to the gradient of solute at the disk times the particle perimeter or 

c _ 1 /• 27rZWao/l 1\ . 

c,in(e/i?) [r: rj' ^ ^ 

where the concentration = co(l + T/Rc) at the cutoff distance is set so that dR/dt = 
for R = Rc, Dc is the solute diffusivity in the matrix and Qq is a microscopic length that 
defines the width of the particle's collection area. For the 2D/2D case, self-similarity with 
one scaling length is consistent with Eq. (|^) since ^/R does not change, and hence the ratio 
of dR/dt for any two particles as given by Eq. is independent of scale as required by Eq. 
(^). The result {R) oc t^^^ then follows from Eq. (|^), self-similarity and Eq. (^. 

In the 3D/2D case, self similarity with one scaling length would require ^ to scale as R 
in Eq. (^). This was assumed by Chakraverty in a mean field treatment with the result 
{R) oc t^/^ which again follows from Eqs. @ and @, and self-similarity. If, however, ^ scales 
as dav, the inter-particle separation, then it must scale in accord with Eq. (|^) and therefore 
^/i? is no longer independent of scale. It would then follow from Eq. (^) that Eq. (^) cannot 
hold, and the particle distribution function cannot be of the form (|^). We present evidence 
in Section |IV| that this is the case, independent of the mean field approximation, and hence 
that self-similarity cannot hold strictly. We point out, however, that the argument of the 
logarithm is expected to change slowly in time, and hence the logarithmic factor itself will 
be changing slowly. Hence, approximate scale invariance could be expected over relatively 
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long periods of time. 



III. THREE DIMENSIONAL PARTICLES ON A ONE DIMENSIONAL MATRIX 

We present in this section the numerical solution of a model system comprising spherical 
three dimensional particles arranged on a closed loop (i.e., a line with periodic boundary 
conditions at the ends), such that each particle can only exchange mass with its two nearest 
neighbors. Such a configuration is intended to model diffusion controlled coarsening of 
three dimensional particles when transport through the matrix takes place preferentially 
along a line (e.g., along a dislocation line). The model, originally introduced by Hunderi 



et al. [Q, is a one dimensional version of the so-called bubble models of grain growth in 
polycrystalline materials. A mean field solution has shown that coarsening proceeds in a 
self-similar fashion, and the growth law for the average particle size has been calculated |T3 



The one dimensional model has been extended to study the existence of self-similarity when 



multiple grain orientations and grain boundary anisotropies are allowed |l6l . 

We consider a set of spherical particles of radii Ri,i = 1, . . . , N, forming a linear chain 
with periodic boundary conditions. Particles are arranged on an evenly spaced grid, such 
that the initial inter-particle separation is one. The rate of change of each particle radius is 
given by 

where M is a mobility coefficient that sets the time scale appropriate for the microscopic 
mechanism responsible for diffusion, and dj.i+i and di^-i are the distances between particles 
i and its two nearest neighbors i + 1 and i — 1 respectively. As the system evolves, some 



particles shrink to zero radius and are removed. Therefore the total number of particles N 
decreases with time whereas both average particle size and average inter-particle distance 
increase. It immediately follows from Eq. (p!0|) that the total volume of the set of particles, 
V = (47r/3) Y.f=i Ri is independent of time. 

We initially place a large number of particles {N = 3 x 10®) on a line, and impose 
periodic boundary conditions RN+i(t) = Ri(t), and dN,i(t) = (ii,Ar(t). The Euler method is 
used to integrate the system of equations ( |T0|) with M = 1 and a step size At = 5. The 
initial condition is a set of randomly chosen radii, uniformly distributed between Rmin and 
Rmax- The value of Rmin is chosen so that the algorithm is stable, and that no particle 
with R > Rmin can shrink to zero radius in At. We have further chosen Rmax = 10 
arbitrarily. The numerical solution proceeds as follows. Given a configuration at time t, 
{Ri(t)}, i = 1, . . . , N(t), Eq. (|T0|) is iterated once for each particle to yield {Ri(t + At)}. Any 
particle for which R{t + At) < Rmin is eliminated, so that only N{t + At) particles remain. 
Links are then redefined so that each particle is connected to its two nearest neighbors but 
preserving their original relative distances. We have studied 420,000 iterations, at the end 
of which 99053 particles were left. We check the accuracy of the integration by monitoring 
the total volume of the set of particles. At very early times, a large number of particles are 
lost and the volume after the first 10,000 iterations decreased by 2.5 %. From then on to 
the end of the calculation, the total volume only changed by 0.3 %. 

Figure |1| shows the scaled distribution of particle radii at three different times. The earlier 
time is slightly before entering the self-similar regime, the other two are arbitrary times 
within it. All other later times agree with these two within the size of the symbols, indicating 
the existence of a self-similar regime. In addition, {R)^ / < d >~ 0.0033 {{d) = dav) changes 

10 



by less that 0.2 % between t = 10^ and the end of the calculation, in agreement with the 
two-length scaling presented in Section |I| (Eq. (Q)). Figure]^ further shows the distribution 
of inter-particle distances scaled with the average particle radius to the third power, also 
consistent with the predictions of Section |I| with regard to the existence of two scaling 
lengths. 

In summary, even though this model has two independent length scales, namely the 
characteristic particle size and inter-particle distance, there is an attractor for the evolution 
in which each distribution function is scale invariant, and each length scale satisfies a well 
defined relationship in the asymptotic limit of long times (although not proportional to each 
other). Finally, Fig. | shows our best estimate of the exponent n, which is very close to the 
theoretical prediction of n = 1/7. 

IV. THREE DIMENSIONAL PARTICLES ON A TWO DIMENSIONAL MATRIX 

Self-similarity and the associated growth exponent are considerably more difficult to 
investigate numerically in this case. Previous numerical approaches in three dimensions have 
considered configurations comprising a large number of point particles and computed growth 
rates by direct summation of the Green's function of the Laplace operator [p!7|-[T9|. We have 
not used this method because the calculation of the Green's function for given boundary 
conditions requires infinite sums over the appropriate images, a procedure that does not 
converge in two dimensions. Furthermore, possible deviations from self-similarity depend 
logarithmic on time, and hence are very weak. As a consequence, numerical integration 
of some particular model would require long time spans to unambiguously discriminate 
between such a dependence and power law growth with some effective exponent. We have 
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therefore focused on the ratio test, namely on whether the ratio of expected growth rates of 
particles of different sizes is independent of scale. If the ratio is not independent of scale, 
then the distribution of particle sizes cannot be scale invariant (i.e., Eq. cannot hold). 
The generality of this procedure is limited by the specific particle distributions that we 
use to apply the ratio test. These distributions are not obtained self-consistently by direct 
solution of a coarsening system. It is therefore conceivable that the failure we find of the 
ratio test for our distributions might not hold for a coarsening system. However, all our 
numerical evidence clearly points to logarithmic factors in the law of boundary motion for 
a two dimensional matrix, and we think it is very unlikely that these factors would cancel 
for a special class of configurations. 

We discuss two cases. First, we present a numerical solution to Laplace's equation in 
a two dimensional square domain with circular disks of fixed radius and concentration at 
the corners and periodic boundary conditions; the configuration is equivalent to an infinite 
square lattice with alternating small and big circles at the lattice sites. The solution es- 
tablishes that, for this simplified configuration of only two types of interacting particles, 
concentration fields in the domain do have a logarithmic factor involving ratios of particle 
sizes to inter-particle distances. Second, we extend this solution to an ensemble of small 
discs with a far field boundary condition of mean-field type, which is also solved for self- 
consistently. The same conclusion holds for this analysis. 

The results will be analyzed in terms of either (^R\Rj directly, or in terms of integrated 
fluxes to particles which are easier to determine numerically. For the latter purpose we define 
Q to be the integrated volume flux transferred from all shrinking particles to all growing 
particles, and q = Q/N; other related quantities are defined in section |1V Ij| . By definition. 
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Q = SopJ^R, t)dR (see Eq. (H) for the definition of Sdj,)- Given the definition of J given 
in Appendix A (Eqs. ( |A9| ) and ( |A3| ) ) , and the result (^, one has, in the self-similar regime, 

g oc {Rf^-^ (11) 

q is constant in the standard 3D/3D case, and inversely proportional to the average particle 
radius for 2D/2D. 

In general, g is a function of the radius and location of the centers of all the particles, 
and of the position of the center of the outer boundary xq (assumed, for example, spherical) 
and its radius Rq: q = q{Ro, Rf, Xq, Xj). The function g is a homogeneous function of degree 
-1 for a two dimensional substrate allowing for capillarity. This follows from the observation 
that if c(r) is a solution of Laplace's equation for the original configuration satisfying all 
boundary conditions, then (1/A)c(r/A) is a solution for the configuration scaled up uniformly 
by A. Hence all gradients scale as 1/A^, and integrated fluxes to each particle and therefore 
g by 1/A. Now, starting with a given spatial configuration, consider scaling up all particle 
radii by a factor s, and all centers of particle positions by a different factor t. Then, 

q{sRutXj) = q (^sRu s^^j^ = ^9 (^Ri', ^ ^j^ = -^f(t/s), (12) 

where /(I) = q{Ri; xj), and, for simplicity, we have taken the origin of the outer boundary 
to lie at xo = 0. The function / depends on the ratio of s and t only. As a consequence, it is 
sufficient to consider rescaling particle distances at fixed particle radii (s = 1) to obtain the 
scaling of q, and, if self similarity holds, to obtain the growth law for the average particle 
size through Eq. (^). 
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A. Two interacting particles in a Laplacian field 

The first case considered is shown in Fig. |. Two discs of radii Ra and and concentra- 
tion Ca and Cb, are placed on the corners of a square lattice of side L with periodic boundary 
conditions. We have considered the case = = 3 and Cf, = —l,Rb = 2 with values 

of L ranging from L = 50 to L = 500. In order to model quasi-static diffusive transport in 
the matrix Laplace's equation is solved in the interior region subject to the boundary con- 
ditions specified at the discs boundaries and periodic boundary conditions otherwise. The 
computational domain is discretized in A''^ evenly spaced elements in each direction. The 
values of Np are adjusted for each L in order to obtain a reasonably accurate solution, and 
range from Np = 500 to Np = 2500. Laplace's equation for the concentration field has been 



solved with a Successive Over-relaxation (SOR) method with Chebyshev acceleration ||20 

Figure || shows the dependence of q with L at constant radii and concentration of the 
particles obtained from the direct solution to Laplace's equation in the square domain. The 
fit indicates the presence of a logarithmic factor in q in the solution, in agreement with the 
results of the mean field calculation presented in Section ^ It shows, for this simple case, 
that the function / in Eq. ([T2|) has an approximate logarithmic dependence on its argument 
t/s. When combined with Eq. ([Til), the result indicates that the coarsening exponent cannot 
be strictly 1/4. Furthermore, since L is proportional to dav, the inter-particle spacing, the 
results indicate that ^ in the mean field expression (Eq. (P)) is not proportional to R but 
rather (approximately) to dav If this is so, it follows that the ratio test cannot be satisfied, 
the reduced particle size is not scale invariant and strict self-similarity does not hold. 
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B. Interacting particles in a Laplacian field 

To investigate the possibility that the growth rates of individual particles all contain the 
same logarithmic factor which therefore cancels out when ratios are formed and allows the 
ratio test to be satisfied, we have extended the calculations to a fairly large collection of 
small interacting particles. The results reported below again indicate that the ratio test fails 
and therefore that strict self-similarity does not hold. 

For numerical reasons, we have chosen elemental square particles of random side ranging 
from the grid spacing h to 5h. This eliminates the need to resolve circular contours, but 
introduces additional anisotropics (and singularities) into the solution. However, we have 
verified that both effects are negligible in the range of parameters used in our calculations. 
Furthermore, c(r) becomes spatially isotropic at distances ~ 3h from the square. 

A large collection of such squares {N = 50) has been placed at random within the inner 
fourth of a square lattice of side L. We have randomly assigned the value Cp = +1 and 
Cp = —1 with equal probability (see Fig. Laplace's equation for the concentration is 
solved in the outer region with the boundary condition c = Coo (uniform) on the outer 
boundary. The total flux through the outer boundary Jqo = — J n ■ Vcdl is then computed. 
For an arbitrary choice of Coo, ^oo 7^ 0. An iterative procedure is then performed by adjusting 
the value of Cqo until Jqo = 0. Such a procedure is intended to model a system of interacting 
precipitate particles embedded in a Laplacian field, chosen, as usual, so that the total mass 
is conserved. 

Once a self-consistent solution for a given L has been found, the linear dimensions of the 
system are scaled up while keeping h and Cp constants and the entire procedure is repeated. 
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Figure ^ shows the value of Coo as a function of L for the configuration shown in Fig. |^. There 
is an evident dependence of on InL. The three parameter fit is needed to account for 
a non-zero value of the average concentration. We have also performed a similar numerical 
computation involving discs of the same size and fixed concentration (so that the average 
concentration can be set to zero). A two parameter fit to the resulting function Cqo vs. L is 
equally good. 

We have repeated all our calculations for the related case that involves a large system 
with periodic boundary conditions, with the elemental squares being uniformly distributed 
throughout the entire computational domain. Identical results have been obtained, the 
details of which will not be presented here. 

We apply the ratio test in the following form: Eq. (j^) implies that if a and j3 are any two 
sets of particles each with specified radii, then the ratio {Qa) / (Q/?) is independent of scale, 
where {Qa) is the expected value of the integrated flux to particles in set a, and similarly 
for {Qp). Thus 

(Q„) = 2vri?- {Ri\R^) = 271 {Rf ^ ^'Mx^) (13) 

and similarly for {Qp) so that in the ratio, the prefactors to the sum cancel and the sums 
are independent of scale {xi here is the reduced particle radius). 

Figures ^ and ^ show the results of our investigation for the configuration of Fig. For 
these cases, it is clear that the ratio test fails. In view of these results, it seems very unlikely 
that there is a special class of configurations for which the ratio of the expected values of the 
Q's for any two subset of particles (each with a specified set of i?j's) could be independent 
of scale. If this is so, then the reduced particle size distribution in the 3D/2D system is not 
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independent of scale and self-similarity does not hold. 
We summarize our principle results: 

1. When the dimensionality of the precipitate particles and of the diffusion matrix differ, 
self-similarity with one scaling length is not possible as it violates mass conservation. 
Scaling with two or more lengths is possible. 

2. In the case of three dimensional particles embedded in a one dimensional diffusion 
matrix, we present evidence of scaling with two lengths; the average particle size {R), 
and the inter-particle spacing {d). The growth laws are (R) oc t^^'^, (d) oc t^/^, and the 
distributions of the corresponding reduced variables become independent of time. 

3. If the distribution of reduced particle size x — R/ (R) is independent of scale, we 
have shown that the expected growth rate of particles of a given radius factors into a 
product of a function of time only (i.e., d (R) / dt) and a function of x only. It follows 
that the ratio of the expected growth rates of particles of two different sizes, or of two 
groups of particles of different sizes, is independent of scale. This ratio test becomes 
a necessary condition for self-similarity. 

4. We have presented evidence that coarsening of three dimensional particles connected 
by a two dimensional diffusion matrix is not self-similar. In particular, we have pre- 
sented evidence from both mean field arguments and from direct numerical solution of 
particles in a Laplacian field, that the ratio test fails. The failure is weak (logarithmic) 
and hence over short periods coarsening may appear self-similar with {R) cx t^/^. 
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APPENDIX A: EXPECTED RATE OF CHANGE OF PARTICLE RADIUS 



We start from the continuity equation for the particle radius distribution 

where n{R,t)dR is the number of particles per unit "area" in the system and J = 
n{R,t)u{R,t) in which u is the average velocity of particles along the R axis, that is 
u = (^R\R^, or the expected value of dR/dt given the value of R. 
Now assume that the particle distribution scales so that 

niR,t) = ^P{x), (A2) 

where P{x) with x = R/ {R) is normalized on x; thus expressed is a product of functions of 
t only and x only. We seek the functional forms of J and u. 

To this end, we first convert Eq. ([All) to express it in terms of x and t. We have. 



^dn\ / dn\ / dn\ / dx\ 

dn\ X d (R) I dn 
dt) ~ {R) dt \dx 



-Tlv^fl^l (A3) 



and 

( dJ\ ( dJ\ ( dx\ 1 ( dJ^ 



dRi. \dRl. \dRl. {R) \dx 



(A4) 



Substituting these expressions into (|A1|) we find 



as the equation of continuity in terms of x and t. Now substitute ( |A3| ) into the right hand 
side of ( [A5D and then add and subtract the term {d{R) /dt){N/ {R))P which allows the 
derivatives to be combined with the result, 
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'dJ\ _d{R) N d{xP) pdN 
^dx j , dt {R) dx dt 



(A6) 



But (r^A = aN {R)^^ by virtue of where a is a constant. Therefore, 



1 dN D„d (R) ^ 
+ .1 \. = 0. 



N dt {R) dt 



Using this relation allows terms of ( [A 61 ) to be combined to yield 



(A7) 



1 d{xP) 
Dp dx 



+ P 



(A8) 



which shows that the left hand side is a product of functions of time and x only. This 
expression may be integrated with the added arbitrary function of time chose to cause J to 



vanish at infinity. The result is, 



J = nu = —N 



X 



-P 



P{x')dx' 



(A9) 



which shows that J and hence u are products of functions of t only and of x only. 
One can go farther by solving Eq. (|A9|) for u and using (|A3|) again to get 



u{x, t) 



R\R 



d{R) 
dt 



G{x) 



(AlO) 



where 



G{x) = x--^ P{x')dx'. 

P J X 



(All) 



Therefore the decomposition (^) or ( [A10| ) is a necessary condition for the scaling assumption 
of the particle size distribution. 
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FIGURES 




R/<R> 



FIG. 1. Scaled distribution of particle radii at three different times as indicated, for three 
dimensional spherical particles on a line. At the earliest time the system is still approaching a 
scale invariant state. 
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400 




FIG. 2. Distribution of inter-particle separations, P{d), scaled by the average particle size (R)^ 
for the times shown. 
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FIG. 3. Logarithmic time derivative of the average particle radius (R) (t) plotted versus 1/ (R). 
For power law growth, the logarithmic derivative asymptotes to a constant equal to the value of 



the exponent n. The arrow on the left axis indicates the exact position of the value 1/7. 
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k 1 

FIG. 4. Two dimensional configuration used in the solution of V^c = 0; the side of the compu- 
tational cell is L = 100. Two discs of radius Ra = S (left bottom) and Rb = 2 (right bottom) are 
placed at the corners of the square domain with periodic boundary conditions. The concentration 
field satisfies Laplace's equation outside the discs, and is constant inside and equal to = 1 and 
Ch = —1. We show in grey scale the iso-concentration lines. 
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FIG. 5. Average rate of volume transfer from particle b to particle a as a function of the system 
size L. The remaining parameters have been kept constant and equal Ra = 3, i?b = 2, Ca = 1 
and Cb = —1. The numerical solution is indicated by the circles, and the solid line is a fit to the 



logarithmic function shown. 
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FIG. 6. Example of the configuration used in the numerical solution of V^c = subject to 
self-consistent mean field boundary conditions at the outer boundary. Shown are 50 elemental 
domains of random size (uniformly distributed between h and 5h, with h the lattice spacing). The 
composition of each domain is chosen randomly as ±1 with equal probability. The 50 small elements 
have been placed at random (uniformly distributed) within the inner quarter of the computational 
domain. We fix Cqo at the outer boundary {x = 0,x = L,y = and y = L) and compute the 
total flux through this boundary. As described in the text, Cqo is then adjusted until the total flux 
vanishes. The solution found is shown in the figure in grey scale. 
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FIG. 7. Mean field concentration Coo as a function of system size L for the configuration shown 
Fig. |. 
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FIG. 8. Ratio test for the configuration shown in Fig. ^. Qa is the average flux to particles 
of size ah (the average taken over the configuration). As shown by the figure, the ratio is not 
independent of L, but includes a logarithmic factor as discussed in the text. 
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FIG. 9. Ratio test for the configuration shown in Fig. ^. Qa is the average flux to particles of 
size ah (the average taken over the configuration). 
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